The embodiments of the invention in which an exclusive property or privilege is claimed are 
defined as follows: 

1. A machine-processing method for computing a property of a mathematically modelled 
physical system, the steps comprising: 

a) reading, via a machine processing unit, input data including a value for each identified 
ordered coefficient of a first polynomial p(x) representing said property, said polynomial p(x) 
being expressed as p(x) = J E i (Pj -*0 where/ = 0 to fl, a value of a quantity x, a value of a 
predetermined x h and a value of a predetermined p(x) correspondingly paired with said 
predetermined x& 

b) building, via said machine processing unit, a value of a second polynomial c(x) 
having ordered coefficients, said second polynomial c(x) being expressible as: 
c(x) = £(Cfc • x k ) where k = 0 to (n - 1) so that said first polynomial p(x) is expressible as: 
P(x)=p(xi) + {x-xi}< comprising the steps of: 

i) determining, via said machine processing unit, a value for each ordered 
coefficient of said second polynomial c(x) by generating a plurality of machine processing unit 
signals to determine each said ordered coefficient of said second polynomial c(x) from: 
Ck = H(Pm+j) • *() where j = 0 to (n - 1 - k) ; 

ii) determining, via said machine processing unit, a value of said second 
polynomial c(x) by generating a plurality of machine processing unit signals to determine: 
c(x) - £(C* • x*) where it = 0 to (n - 1) ; 
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c) constructing, via said machine processing unit, a value of said first polynomial p(x) by 
generating a plurality of machine processing unit signals to determine: 

P(x) = P(*i) + i x - x i } • c ( x ) ; and 

d) outputting, via said machine-processing unit, said value of the first polynomial p(x) 
5 representing said property of the mathematically modelled physical system. 



2. The machine-implementable method of claim 1, wherein a difference between x and x, is 
sufficiently small to achieve a desired accuracy of a final computation of said machine 
representation of a numerical value of said first polynomial p(x). 



'^10 

jc 



3. The machine-implementable method of claim 2 wherein the step of reading said input 
data comprises reading, via said machine processing unit, said input data from a 
machine-readable medium. 



15 4. The machine-implementable method of claim 3 wherein said ordered coefficients of said 
second polynomial c(x) are computed from a mathematical expression derivable from: 

C k = Z(iVi+fl * */) where 7 = 0 to (« - 1 - 



5. The machine-implementable method of claim 4 wherein said mathematical expression is 
20 a mathematical recurrence expression. 
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6. The machine-implementable method of claim 5 wherein said mathematical recurrence 
expression is a forward mathematical recurrence expression. 

7. The machine-implementable method of claim 5 wherein said mathematical recurrence 
expression is a backward mathematical recurrence expression. 

8. The machine-implementable method of claim 7 further adapted to implement said 
backward mathematical recurrence expression by comprising further steps for: 

e) equating, via said machine-processing unit, a value of a highest ordered coefficient of 
said second polynomial c(x) to a value of an identified highest ordered coefficient of said first 
polynomial p(x) by generating a plurality of machine processing unit signals to determine: 
C n -i = Pn\ and 

f) computing, via a machine processing unit, a value for each lower ordered coefficient 
of said second polynomial c(x) by generating a plurality of machine processing unit signals to 
determine: C k -\ = x t -C k + P k where fc=(w- 1) to 1 . 

9. The machine-implementable method of claim 8 wherein said predetermined x £ is selected 
from a set of predetermined values of x h 
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10. The machine-implementable method of claim 9 wherein said predetermined x t is a 
closest member of said set to said identified x. 

11. The machine-implementable method of claim 10 wherein said step of determining a 
5 value of said second polynomial c(x) is computed by using Horner's Rule. 

12. The machine-implementable method of claim 11 for determining a value of a 
O denominator polynomial q(x) having identified ordered denominator coefficients, said 
g denominator polynomial q(x) being expressible as: q(x) = 2(2/ • *0 where/ = 0 to ro, 
\jlo comprising further steps of: 

I, g) adapting said input data to further include a value for each identified ordered 

fh denominator coefficient of said denominator polynomial q(x) 9 a value of a predetermined qfx) 

\& 

Q correspondingly paired with said predetermined x& and 

h) determining, via said machine processing unit, a value of another polynomial dfx) 
15 having ordered denominator coefficients, said another polynomial dfx) being expressible as: 
d(x) = 2(Da; -x*) where fc = 0 to (m - l)so that said denominator polynomial #(x) is expressible 
as: g(x) = #(x;) + {x -x/} • d(x) 5 and a value for the said denominator polynomial is resolved. 

13. The machine-implementable method of claim 12 wherein the first polynomial pfx) is a 
20 numerator polynomial pfx), and p(x)-p(x) is computed, and pfx) is not read. 
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14. The machine-implementable method of claim 13 for determining a value of a rational 
function r(x) being expressible as a quotient of said numerator polynomial p(x) and said 

p(x) 

denominator polynomial q(x) expressed as r(x)= , comprising further steps of: 

5 i) adapting said input data to further including a value of a predetermined K x d 

correspondingly paired with said predetermined x y ; and 

j) constructing, via said machine processing unit, a value of said rational function r(x) by 
generating a plurality of machine processing unit signals to determine: 

r(x) = r(x ,) . (1 - ?(x) ) + ^ . 



10 



15. The machine-implementable method of claim 14 wherein said rational function r(x) is an 



Q approximation to a Bessel function. 



16. The machine-implementable method of claim 14 wherein said rational function r(x) is an 
15 approximation to an error function (ERF). 



17. The machine-implementable method of claim 14 wherein said rational function r(x) is an 
approximation to a complementary error function (ERFC). 
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18. The machine-implementable method of claim 14 wherein said rational function r(x) is an 
approximation to a log gamma function (LGAMMA). 



19. The machine-implementable method of claim 11 or 14 wherein all values are machine 
5 representations of some numerical value, said machine processing unit is a computer processing 
unit, each machine representation is a binary representation operable with said computer 
processing unit, and machine-readable medium is a computer-readable medium. 



20. A computer-readable program product having computer executable instructions for 
'iio instructing a computer, said instructions tangibly embodying the machine-implementable 
method of claim 19. 



21. A computer-readable mathematical software routine library including computer 
executable instructions for instructing a computer, said instructions tangibly embodying the 
15 machine-implementable method of claim 19. 



22. The computer-readable mathematical software routine library of claim 21 wherein said 
library is operatively associated with a software programming language. 



CA920000040US1 32 



23. A machine for computing a property of a mathematically modelled physical system, the 
steps comprising: 

a) reading, via a machine processing unit, input data including a value for each identified 
ordered coefficient of a first polynomial p(x) representing said property, said polynomial p(x) 
being expressed as p(x) = • x j ) where j = 0 to n 9 a value of a quantity x, a value of a 
predetermined x /? and a value of a predetermined p(x) correspondingly paired with said 
predetermined x { \ 

b) building, via said machine processing unit, a value of a second polynomial c(x) 
having ordered coefficients, said second polynomial c(x) being expressible as: 
c{x) = Y*{Ck • x*) where A: = 0 to (n - 1) so that said first polynomial p(x) is expressible as: 
p(x) = p(x;) + {x -x i} - c(x) 9 comprising the steps of: 

i) determining, via said machine processing unit, a value for each ordered 
coefficient of said second polynomial c(x) by generating a plurality of machine processing unit 
signals to determine each said ordered coefficient of said second polynomial c(x) from: 
C k = 1i(P(k+i+f> • */) where j = 0 to (n - 1 - k) ; 

ii) determining, via said machine processing unit, a value of said second 
polynomial c(x) by generating a plurality of machine processing unit signals to determine: 
c(x) = ^(Ck • x k ) where k = 0 to O - 1); 

c) constructing, via said machine processing unit, a value of said first polynomial p(x) by 
generating a plurality of machine processing unit signals to determine: 
p(x) = p(xi) + {x - xi } • c(x) ; and 



CA920000040US1 33 



d) outputting, via said machine-processing unit, said value of the first polynomial p(x) 
representing said property of the mathematically modelled physical system. 

24. The machine of claim 23 wherein a difference between x and x f - is sufficiently small to 
achieve a desired accuracy of a final computation of said machine representation of a numerical 
value of said first polynomial p(x). 

25. The machine of claim 24 wherein said means for reading said input data comprises 
means for reading, via said machine processing unit, said input data from a machine-readable 
medium. 

26. The machine of claim 25 wherein said ordered coefficients of said second polynomial 
c(x) are computed from a mathematical expression derivable from: 

C k = • *{) where j = 0 to (n - 1 - k) . 

27. The machine of claim 26 wherein said mathematical expression is a mathematical 
recurrence expression. 

28. The machine of claim 27 wherein said mathematical recurrence expression is a forward 
mathematical recurrence expression. 
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29. The machine of claim 27 wherein said mathematical recurrence expression is a backward 
mathematical recurrence expression. 

30. The machine of claim 29 further adapted to implement said backward mathematical 
recurrence expression by further comprising: 

e) means for equating, via said machine processing unit, a value of a highest ordered 
coefficient of said second polynomial c(x) to a value of an identified highest ordered coefficient 
of said first polynomial p(x) by generating a plurality of machine processing unit signals to 
determine: C«-i = Pn\ and 

f) means for computing, via said machine processing unit, a value for each lower ordered 
coefficient of said second polynomial c(x) by generating a plurality of machine processing unit 
signals to determine: C*-i = -C k + P k where k = (» - 1) to 1 . 

31. The machine of claim 30 wherein said predetermined jc, is selected from a set of 
predetermined values of x t . 

32. The machine of claim 30 wherein said predetermined x, is a closest member of said set to 
said identified x. 
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33. The machine of claim 32 wherein the determining means for determining a value of said 
second polynomial cfx) is computed by using Horner's Rule. 

34. The machine of claim 33 for determining a value of a denominator polynomial q(x) 
having identified ordered denominator coefficients, said denominator polynomial q(x) being 
expressible as: q(x) = 2(6/ ' *0 where j = 0 to m 9 comprising further steps of: 

g) adapting said input data to further include a value for each identified ordered 
denominator coefficient of said denominator polynomial qfx), and a value of a predetermined 
q(Xj) correspondingly paired with said predetermined x t \ and 

h) determining, via said machine processing unit, a value of another polynomial d(x) 
having ordered denominator coefficients, said another polynomial dfx) being expressible as: 
d(x) = YdiPk • x h ) where k = 0 to (m - 1 )so that said denominator polynomial qfx) is expressible 
as: q(x) = q(xi) + {x -Xi} • 5 and a value for the said denominator polynomial is resolved. 

35. The machine of claim 34 wherein the first polynomial pfx) is a numerator polynomial 
pfx), and p(x)-p(x) is computed, and ^(xj is not read. 

36. The machine of claim 35 for determining a value of a rational function rfx) being 
expressible as a quotient of said numerator polynomial pfx) and said denominator polynomial 

p(x) 

qfx) expressed as rfx)= , comprising further steps of: 
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i) adapting said input data to further including a value of a predetermined K x d 
correspondingly paired with said predetermined x t ; and 

j) constructing, via said machine processing unit, a value of said rational function r(x) by 
generating a plurality of machine processing unit signals to determine: 

fix) = r{ Xi ) . (1 - ^ ) + ?(x) . 

37. The machine of claim 36 wherein said rational function is an approximation to a Bessel 
function. 

38. The machine of claim 36 wherein said rational function is an approximation to an error 
function (ERF). 

39. The machine of claim 36 wherein said rational function is an approximation to a 
complementary error function (ERFC). 

40. The machine of claim 36 wherein said rational function is an approximation to a log 
gamma function (LGAMMA). 
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41. The machine of claim 33 or 36 wherein all values are machine representations of some 
numerical value, said machine processing unit is a computer processing unit, each machine 
representation is a binary representation operable with said computer processing unit, and said 
machine-readable medium is a computer-readable medium. 

42. A machine having a computer-readable program product having computer executable 
instructions for instructing a computer to embody the machine of claim 41 . 

43. A machine having a computer-readable mathematical software routine library including 
computer executable instructions for instructing a computer to embody the machine of claim 41 . 

44. A machine having the computer-readable mathematical software routine library of claim 
43 wherein said library is operatively associated with a software programming language. 
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